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Abstract: We analyze the time evolution of a one-dimensional quantum sys- 
tem with an attractive delta function potential whose strength is subjected to a 
time periodic (zero mean) parametric variation r/(t). We show that for generic 
Tj(t), which includes the sum of any finite number of harmonics, the system, 
started in a bound state will get fully ionized as t — > oo. This is irrespective of 
the magnitude or frequency (resonant or not) of rj(t). There are however excep- 
tional, very non-generic rj{t), that do not lead to full ionization, which include 
rather simple explicit periodic functions. For these rj{t) the system evolves to 
a nontrivial localized stationary state which is related to eigenfunctions of the 
Floquet operator. 

Key words. Ionization; Resonances; Time-dependent Schrodinger equation 



1. Introduction and results 

We are interested in the qualitative long time behavior of a quantum system 
evolving under a time dependent Hamiltonian H(t) = Hq + Hi(t), i.e. in the 
nature of the solutions of the Schrodinger equation 



ihd t i) = [H Q + H x {t)]i } (i) 

Here ip is the wavefunction of the system, belonging to some Hilbert space H, Ho 
and Hi are Hermitian operators and equation(l) is to be solved subject to some 
initial condition ipQ. Such questions about the solutions of ([!]) belong to what 
Simon M calls "second level foundation" problems of quantum mechanics. They 
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are of particular practical interest for the ionization of atoms and/or dissociation 
of molecules, in the case when Hq has both a discrete and a continuous spectrum 
corresponding respectively to spatially localized (bound) and scattering (free) 
states in R d . Starting at time zero with the system in a bound state and then 
"switching on" at t = an external potential H\(t), we want to know the 
"probability of survival", P(t), of the bound states, at times t > 0: P(t) = 
J2j KVK*)' u i)l 2 ' where the sum is over all the bound states Uj ||,||,[|,||,|6|,||,[|. 

This problem has been investigated both analytically and numerically for the 
case H\(t) = r)(t)Vi(x) with ij(t) = rs'm(ujt + 9) and V\ a time independent 
potential, x € R d . When ui is sufficiently large for "one photon" ionization to 
take place, i.e., when huj > —Eq, Eq the energy of the bound (e.g. ground) state 
of Hq and r is "small enough" for Hi to be treated as a perturbation of Hq then 
this is a problem discussed extensively in the literature (&0j). Starting with the 
system in its ground state the long time behavior of P(t) is there asserted to 
be given by the P(t) ~ exp[— rpt]. The rate constant rp is computed from first 
order perturbation theory according to Fermi's golden rule. It is proportional to 
the square of the matrix element between the bound and free states, multiplied 
by the appropriate density of continuum states in the vicinity of the final state 
which will have energy fko — E a mMmu^. 



Going from perturbation theory to an exponential decay involves heuristics 
based on deep physical insights requiring assumptions which seem very hard to 
prove. It is therefore very gratifying that many features of this scenario have 
been recently made mathematically rigorous by Soffer and Weinstein || (their 
analysis was generalized by Soffer and Costin J5J). They considered the case 
when Hq = — V 2 + Vq(x), x £ R 3 , Vq compactly supported and such that there 
is exactly one bound state with energy — ui (from now on we use units in which 
h = 2m =1) and a continuum of quasi-energy states with energies k 2 for all 
k G R 3 . The perturbing potential is Hi(t) — r cos(ujt)Vi(x) with V%(x) also of 
compact support and satisfying some technical conditions. They then showed 
that for lo > loq and r small enough there is indeed an intermediate time regime 
where P(t) has a dominant exponential form with the Fermi exponent Pp. This 
regime is followed for longer times by an inverse power law decay. Some of these 
restrictions can presumably be relaxed but the requirement that r be small is 
crucial to their method which is essentially perturbative. 

The behavior of P{t) becomes much more difficult to analyze when the 
strength of H\{t) is not small and perturbation theory is no longer a useful 
guide. This became clear in the seventies with the beautiful experiments by 
Bayfield and Koch, c.f. for a review, on the ionization of highly excited 
Rydberg (e.g. hydrogen atoms) by intense microwave electric fields. These ex- 
periments, showed quite unexpected nonlinear behavior of P(t) as a function of 
the initial state, field strength E and the frequency uj. These results as well as 
other multiphoton ionizations of hydrogen atoms have been (and continue to be) 
analyzed by various authors using a variety of methods. Prominent among these 
are semi-classical phase-space analysis, numerical integration of the Schrodinger 
equation, Floquet theory, complex dilation, etc. While the results obtained so far 
are not rigorous, they do give physical insights and quite good agreement with 
experiments although many questions still remain open even on the physical 
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In addition to the above experiments on Rydberg atoms there are also many 
experiments which use strong laser fields to produce multiphoton (u < —E$) ion- 



ization of multielectron atoms and/or dissociation of molecules [116 , 17 . These 



systems are more complex than Rydberg atoms and their analysis is correspond- 
ingly less developed. One unexpected result of certain studies is that an increase 
in the intensity of the field may reduce the degree of ionization, i.e., P(t) can be 
non monotone in the field strength E at large values of E. This phenomenon, 
which is often called "stabilization" , can be observed in some numerical simu- 
lations, analyzed rigorously in some models and is claimed to have been seen 
experimentally c.f. g and E|,^^,^|. 

It turns out that many features observed for Rydberg atoms and also stabi- 
lization are already present in a simple model system which we have recently 
begun to investigate analytically [p2[|2^ ,p^|. The results described there show 
that phenomenon of ionization by periodic fields is very complex indeed once 
one goes beyond the perturbative regime. This will become even clearer from 
the new results about this model presented here. 



2. The model 

We consider a very simple quantum system where we can analyze rigorously 
many of the phenomena expected to occur in more realistic systems described by 
(Q). This is a one dimensional system with an attractive delta function potential. 
The unperturbed Hamiltonian Hq has, in suitable units, the form 



Ho 



-—— — 2S(x), — oo < x < oo. 
ax 1 



(2) 



The zero range (delta-function) attractive potential is much used in the litera- 
ture to model short range attractive potentials [p5 , 26 , 2^,p8|| . It belongs, in one 
dimension, to the class K\ Hq has a single bound state Ub(x) — e _ ' x ' with 
energy — luq — — 1. It also has continuous uniform spectrum on the positive real 
line, with generalized eigenfunctions 



u{k, x) 

and energies k 2 . 

Beginning at t = 
we have 



1 



/2vr 



,/7.'./- 



1 



J\kx\ 



l + t|fc| 



-oo < k < oo 



0, we apply a parametric perturbing potential, i.e. for t > 



H(t) = H -2 V (t)S(x) (3) 

and solve the time dependent Schrodinger equation (0) for ip(x, t), with tp(x, 0) = 
ij)o(x). Expanding ip in eigenstates of Hq we write 

ip(x,t) = 6»(t)u fc (a;)e it 

9(k 1 t)u(k 1 x)e- lk2t dk (t>0) (4) 
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with initial values 9(0) — 9q, 0(k,O) — Oo(k) suitably normalized, 



(rPo^a) = \9 q \ 2 + / \0 o (k)\ 2 dk = l (5) 



We then have that the survival probability of the bound state is P(t) = 
\9(t)\ 2 , while \0(k, t)\ 2 dk gives the "fraction of ejected particles" with (quasi-) 
momentum in the interval dk. 

This problem can be reduced to the solution of an integral equation in a single 
variable H,H. Setting 

Y(t) =ip(x = 0,t)r}{t)e u (6) 

we have 

9{t) = 9 Q + 2i J* Y(s)ds , (7) 
0(k,t) =O (k) + 2|fc|/[v^F(l-i|A:|)] f* Y(s)e^ 1+k2 >ds . (8) 
Y(t) satisfies the integral equation 

Y(t) = r](t) jl(t) + J [2i + M(t - t')]Y{t')dt'\ = v(t)(l(t) + (2i + M) * f) 

(9) 



where the inhomogeneous term is 



and 



7T 

with 



m =9 + ' r M^Be-^)'*, 

V/27T in l+ik 



2i f°° M 2 e - ls (i+« 2 ) 1 + i f°° e~ iu , 



ii 



f*9= f(s)g(t-s)ds 



o 



In our previous works we considered the case where Oo(k) — and ri(t) is a 
finite sum of harmonics with period 2ituj~ 1 . In particular, we showed in ||23| how 
to compute the survival probability P(t) as a function of the strength r and 
frequency u) when rj(t) = rsvauit. Here we study the general periodic case and 
write 



3=0 

Our assumptions on the Cj are 
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(a) O^r/e L°°(T) 

(b) C = 0_ 

(c) C-i = Cj 

Genericity condition (g). Consider the right shift operator T on ^(N) given 

by 

T(Ci, C2, C n , ...) = (C2, C3, C n+ i, ...) 

We say that C G ^(N) is generic with respect to T if the Hilbert space 
generated by all the translates of C contains the vector e\ = (1, 0, 0..., ) (which 
is the kernel of T): 



ei e \j T n C (10) 

(where the right side of ( [To| ) denotes the closure of the space generated by the 
T n C with n > 0.) This condition is generically satisfied, and is obviously weaker 
than the "cyclicity" condition h(N) G V^Lo^™^ = {^}: which is also generic 
p9[ (Appendix B discusses in more detail the rather subtle cyclicity condition). 

An important case, which satisfies (|l0|), (but fails the cyclicity condition) 
corresponds to r\ being a trigonometric polynomial, namely C ^ but C n = 
for all large enough n. (We can in fact replace ei in (|Io| ) by with any k > 1.) 
A simple example which fails ( |l0| ) is 

A — cos(ujt) 

n(t) = 2r\ , v ' r 11 

,W l + A 2 -2Acos(wt) V ; 

for some A € (0,1), for which C n = —rX n for n > 1. In this case the space 
generated by T n C is one-dimensional. We will prove that there are values of r 
and A for which the ionization is incomplete, i.e. 9(t) does not go to zero for 
large t. 



3. Results and Remarks 

Theorem 1 Under assumptions (a). ..(c) and (g), the survival probability P(t) 
of the bound state Ub, \6(t)\ 2 tends to zero as t — > 00. 

Theorem 2 For ijJo( x ) — u b{%) there exist values of X, u> and r in ([/"?[), for 
which \0(t)\ -/-* as t — > 00. 

Remarks. 

1. Theorem [j] can be extended to show that J D \tp(x, t)\ 2 dx — * for any 
compact interval Dei This means that the initially localized particle really 
wanders off to infinity since by unitarity of the evolution J 9 \ip(x,t)\ 2 dx = 1. 
Theorem [2] can be extended to show that for some fixed r and ui in ( |ll| ) there 
are infinitely many A, accumulating at 1, for which Q(t) -/-> 0. In these cases, it 
can also be shown that for large t, 6 approaches a quasiperiodic function. 
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2. While Theorem |l| holds for arbitrary ipo, care has to be taken with the 
initial conditions for Theorem ^. In particular we cannot have an initial state 
such that in @ I(t) = for al t. This would occur, for example, if ipo (x) is 
an odd function of x. In that case the evolution takes place as if the particle 
was entirely free — never feeling the delta function potential. There may also 
be other special ipo for which 9q ^ but for which 9(t) — > as t — > oo. We have 
therefore stated Theorem |^ for the case ipo = tt&. We shall also, for simplicity, 
use this choice of ipo m the proofs of Theorem l[ For this case, which is natural 
from the physical point of view, I(t) = 1 in (3). The extension to general ipo is 
immediate and is given at the end of Section 5. 

3. In we gave a detailed picture of how the decay of 9{t) depends on r and 
oj when r](t) — r sm(cut), 9q = 1. For small r and o> _1 not too close to an integer 
we get an exponential decay with a decay rate r(r, uj) ~ r 2 (i+[w J ) w i icrc [a; -1 J 
is the integer part of uo^ 1 . (For uo > 1, this corresponds to _T ~ JV)- At times 
large compared to -T -1 , \0(t)\ decays as t~ 3 / 2 . The picture becomes much more 
complicated when r is large and/or uj^ 1 is an integer. In particular there is no 
monotonicity in \9(t)\ as a function of r. In [ p4| we proved complete ionization 
for the case where C n = for n > N, N > 1. 

4. We note here that Pillet proved complete ionization for quite general 
Ho under the assumption that Hi(t) is "very random", in fact a Markov process. 
Our results are not only consistent with this but support the expectation that 
generic perturbations will lead to complete ionization for general Ho- This is 
what we expect from entropic considerations — there is just too much phase 
space "out there" . The surprising thing is that even for our simple example one 
can readily find exceptions to the rule. 

We should also mention here the work of Martin et al. J|l],[32| who consider 
the case where Ho has an isolated eigenvalue Eq plus an absolutely continuous 
spectrum in the interval [0,E max \. They show that if the frequency u> of the 
periodic, small, perturbation H\{t) is larger than Eq then the bound state is 
stable. This can be understood in terms of Fermi's golden rule by noting that 
the density of states at the energy Eq + u) > E max is zero so that rp would be 
zero. 

5. There is a direct connection between our results and Floquet theory where, 
for a time-periodic Hamiltonian H(i) with period T — 2tt/uj, one constructs a 
quasienergy operator (QEO) [p||3^,p4|. 

K = -i§- 9+ H(0) 

K acts on functions of x and 9 , periodic in 9, i.e. on the extended Hilbert space 
TL ® £2(5, T~ 1 d9). Let now cp(x, 9) be an eigenfunction satisfying 



K(p = n^, <j>(x, 9 + T) = <p(x, 9) (12) 

then, 

ip(x,t) = e^VOM) 

dt 



is a solution of the Schrodinger equation = H(t)ip. 
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The existence of a real eigenvalue fj, of the QEO with an associated <f)(x, 6) e 
L 2 QSL d (g>S) is thus seen to imply the existence of a solution of the time-dependent 
Schrodinger equation which is, in absolute value, periodic. This shows that, 
for appropriate initial conditions, the particle has a nonvanishing probability 
of staying in a compact domain and thus, for the case considered here, that 
ionization is incomplete. We also note that for each such /i there is actually a 
whole set [i n = [x + nto of eigenvalues of K. 

For the specific model considered here, jl2] ) takes the form 

K(f) = _ a 2 y.0) _ 2(1 + ^e))^ -M = »<t> (is) 

ax z ad 
We can now look for solutions of ( |l3| ) in the form 

with = ±y//j. — nuj. Such a solution is in L 2 only if $t(a n x) < 0, a condition 
which obviously selects different roots X n depending on whether x > or x < 0. 
The requirement that (j)^ be in L 2 (R) leads to a set of matching conditions which 
determine whether such eigenvalues \x can exist. It is easy to see that has to 
be continuous at zero and satisfy the condition 

2^(0", 6) - <^(0+ 6) = 2(1 + 0). 

This implies, after taking the Fourier coefficients of both sides of the above 
equality, the recurrence relation 



y n (2-a+ + a-) = 2^C jyn ^ (14) 

for which a (nontrivial) solution y n € I 2 is sought. This is effectively the same 
equation as d2Q) below which is at the core of our analysis. Complete ionization 
thus corresponds to the absence of a discrete spectrum of the QEO operator 
and conversely stabilization implies the existence of such a discrete spectrum. 
In fact, an extension of Theorem ^ shows that for the initial condition ipo — Uf,, 
ipt approaches such a function with /i = — sq. More details about Floquet theory 
and stability can be found in (3^,Q 

6. We are currently investigating extensions of our results to the case where 
Ho = —V 2 + Vo(x), x £ M. d , has a finite number of bound states and the 
perturbation is of the form rj{t)V\{x) and both Vq and V\ have compact support. 
Preliminary results indicate that, with much labor, we shall be able to generalize 
Theorem^], to generic V\(x). The definition of genericity will, however, depend 
strongly on Vq. 

The physically important case of an external electric dipole field, V\(x) 
I ■ can be transformed into the solution of a Schrodinger equation of the form 
H{t) = —V 2 + Vq(x — g(t)), see This should, in principle, also be amenable 
to our methods but so far we have no results for that case. 

Outline of the technical strategy. The method of proof relies on the 
properties of the Laplace transform of Y, y(p) = CY(p) = J °° e~ pt Y(t)dt. 
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Since the time evolution of -0 is unitary, \0(t)\ < 1. This gives some a priori 
control on Y. For our purposes however it is useful to characterize directly the 
solution of the convolution equation (^). (We restrict ourselves to 6>o(fc) = 
and I(t) = 1 there.) We show that this equation has a unique solution in suit- 
able norms. This solution is Laplace transformable and the Laplace transform y 
satisfies a linear functional equation. 

The solution of the functional equation satisfied by the transform of Y is 
unique in the right half plane provided it satisfies the additional property that 
y(pa + is) is square integrable in s for any po > 0. Any such solution y transforms 
back (by the standard properties of the inverse Laplace transform) into a solu- 
tion of our integral equation with no faster than exponential growth; however 
there is a unique locally integrable solution of this equation, and this solution is 
exponentially bounded. This must thus be our Y. We can thus use the functional 
equation to determine the analytic properties of y(p). 

This is done using (appropriately refined versions of) the Fredholm alterna- 
tive. After some transformations, the functional equation reduces to a linear 
inhomogeneous recurrence equation in I2, involving a compact operator depend- 
ing parametrically on p, see e.g. (|l7|). The dependence is analytic except for a 
finite set of poles and square-root branch-points on the imaginary axis and we 
show that the associated homogeneous equation has no nontrivial solution. We 
then show that the poles in the coefficients do not create poles of y, while the 
branch points are inherited by y. The decay of y(p) when |9(p)| — > 00, and the 
degree of regularity on the imaginary axis give us the needed information about 
the decay of Y{t) for large t. 

4. Behavior of y(p) in the open right half plane H 

Lemma 3 (i) Equation has a unique solution Y £ Lj oc (R + ) and \Y(t)\ < 
Ke Bt for some K, B el. 

(ii) The function y{p) = CY exists and is analytic in Mb = {p ■ > B}. 

(Hi) In Mb, the function y{p) satisfies the functional equation 



y= C^fh + by) (15) 

j=—oo 

with 



(r/) 0) = f(p + tw), hip) = -p- 1 and b(p) = -- (l + y/T^ip) 

The branch of the square root is such that for p € M = {p : 3?(p) > 0}, the 
real part of i/l — ip is nonnegative and the imaginary part nonpositive. 

The straightforward proofs of this lemma are done in Appendix A. (Some of 
the results can also be gotten directly from standard results on the Schrodinger 
operators and on integral equations.) 
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Remark 4 It is clear that the functional equation only links points on 

the one dimensional lattice {p + iLlS\. It is convenient to take po such that 
p = Po + inuj with ^R(po) = 3?(p) and 

9f(po)e[0,«) (16) 

The functions y, h, b in ( [l5|) will now depend parametrically on pq. We set y = 
{yj}j£Z, h = {hj} jeZ , b {bjjjcz with y n = y(p + inu) = y(p) (and similarly 
for h(p) and b(p)). It is convenient to define the operator (Hy) n — b n y n . Let 
(Ty)„ = y n +i be the right shift on ^(Z) (which we denote for simplicity by 1 2) 
and rewrite ( filf ) as 



y= C 3 Tih+ C^Hy=f + Jy (17) 

j— —00 j— — 00 

Proposition 5 For 5R(po) > there exists a unique solution of (iTa) in l-i- This 
solution is analytic in po,^t(po) > 0. Thus y(p) is analytic in p el and inverse 
Laplace transformable there with C^ 1 (y) = Y . 

Proof The proof uses the Fredholm alternative. We first prove the following 
results. 

Lemma 6 The operator J is compact on I2 if Po =/= 0. 

Proof. The proof uses standard compact operator results, see e.g. |3C|l , First 
note that the operator T 3 H is compact. This is straightforward: since bj — > as 
j — ► 00, it follows that H is the norm limit as N — > 00 of the finite rank operators 
defined by (iJjvy)j = bjyj for |j| < N and (HNy)j — otherwise, and thus is 
compact. Since ||T|| = 1, T 3 H is compact too. As the series J2 ( jL^ QO CjT J H 
converges strongly, J is compact. 

Remarks 

1. Note that / S I2 if po 7^ (a straightforward consequence of the fact that 
C and h in (17) are in I2). 

2. The operator J is analytic in pg, except for po = 0, where the coefficients 
have poles, and for an additional value on the imaginary axis (possibly also 0), 
where the coefficients have square root branch points. 

Remark 7 Setting, for po 7^ 0, 

yi = (y/l-i(po + ilu) - l)zi (18) 

the homogeneous equation 

y = Jy (19) 

clearly has a (nontrivial) I2 solution y only if 



00 

(yi - ip + Ioj - lj Zl = - (Ckzi+k + C k zi- k ^ (20) 

fe=i 
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has a (nontrivial) 1% solution z with 

{ - ip +juj - l)^} G l 2 (21) 



Lemma 8 For any r\ under assumptions (a) to (c), ifpo £ H i/iere zs no nonzero 
I2 solution of $2t\ ) such that ^}\ ) holds. 



Proof. To get a contradiction, assume z <E I2, z ^ 0, satisfying (|2l|), is a solution 
of (p0|). Multiplying (|2^) by zj, and summing with respect to I from —00 to +00 
we get 



00 

(Vi-m + iu-i)\z\t 

I— — oo 

00 00 00 00 

= - ^ ^(CkZi+kzJ+CkZ^kzjj = - ^ (c k ziz[^k + CkZi-kTi 



1 — — OC k—l I——00 k—l 

00 00 



I— —00 k—l 



If Po G H the imaginary part of VI — ipo + ^ is negative (see ( |A5| ) and the 
discussion following it) and thus, if some zi is nonzero then the left side of ( p2f ) 
has strictly negative imaginary part, which is impossible since the right side is 
real. 

Proof of Proposition The existence of the analytic solution follows now 
immediately from the analytic Fredholm alternative and the analyticity of the 
coefficients, for po £ EL The fact that {y n } & h together with the stated analyt- 
icity imply that the function £~ 1 y(p) exists and satisfies the integral equation 
of Y, and thus coincides with Y. 



5. Behavior of y(p) in the neighborhood of 3(p) = in the generic case 
Discussion of methods. We start again from relation (|l7|). This has the form 

C • 

y n = iY] — . / . .x V Cjq n+ jy n+j , C = (23) 

^ — ' —inn -V- In. -V- n)in t—' 



-ip + (n + j)uj 



where 



_ [1 + y/1 - ipo + nuj\ 

In — . {^J 

—ipo + noj 

As the imaginary axis $i(pa) — is approached, two types of potential singu- 
larities in the coefficients need attention: the poles in the coefficients due to the 
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presence of and the square root singularities. It will turn out that by can- 
celation effects, the poles play no role, generically. The square root singularities 
will be manifested in the solution y. The st udy of these questions requires further 
regularization of the functional equation (|23| ) . 

It is convenient to separate out the terms in (|2^) which are singular at po = 0. 
Using (from now on) the notation sq = —ipo we have 



.C_„ C_„(l + VT+Ja) , . Cj 
y n = «— ; yo + I " 



s s tr 1 s + (n + 

~ E Cjln+oVn+j, n^O 
C 

y° = 1 E -±- - E cm ( 25 ) 

j^O J j^O 

We break up the proof into two parts, the non-resonant and resonant case. 
We start with the former, 



5.1. The Non-Resonant Case, to 1 $ N. 

Proposition 9 If condition (g) is satisfied, and lo^ 1 $ N, then the solution y 
of ( pq ) is analytic in a small neighborhood of so = 0. 

For the proof we write yo = i/2+Soito, and f° r i^Owe make the substitution 
Vn = Vn + d n Uo, where we will choose d n according to (^6|) in order to eliminate 
u from all equations with fi^O. 

Lemma 10 (i) For s S R there exists a unique solution d G l^ifL \ {0}) of the 
system 

d n = -C-n(l + VI + so) - Yj Ck-nq k d k , n^O (26) 

fc^O 

This solution is analytic at sq = 0. 



(ii) With this choice of d, the system (25) becomes 



so + C i1j d o ) u o = fo - E C iWi ( 27 ) 



where 



fo = ~n + 1 y. — 7~~ ' = lC - n — o + z 2^ — n - ( 28 ) 

2 ^— 'so+jw 2s f-fso + kuj 

j#o 

(mj For small So we have Yljjto ^jljdj 7^ 0, and the system $2y has a unique 
solution with v <E l%ifL \ {0}), and v n , uq are analytic at sq = . 
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Proof, (i) Equation (gff is of the form (I - J')d = c' in 1 2 {Z \ {0}), where 
c' n = -(1 + VTTJ^)C n and 

(J'd) n = -J2 C k-nQkd k , (n?0) 

fe#0 

We show first that Ker(I — J') — {0}. Indeed, let d — Ad and set Dk = qkdk- 
Then we see that 



q n - 1 D n + Y,C k -nD k = Q (29) 

fe^O 

and, by multiplying with D n and summing over n we get 

^ q-^Dn] 2 + C k - n D k D^ = (30) 

Note that, because C_„ = C„, the following quantity is real: 



£ ft^^^ = ^ C„_ fc ^ fc Z?„ = J2 C k - n D k D n (31) 

n,fe^0 n,k^0 n,k=L0 

implying that 
with (cf. (H)) 



g„ = -l + VT+ s + 



Let N = -(1 + so)^ 1 e R. Obviously g^ 1 e R for n > N Q while for n < N 
we have, by Lemma [| (iii) 

^k; 1 ) < o 

Thus it is necessary that D n = for all n < Nq. 

Assume D ^ 0. Let N G N be such that £>„ = for all n < N and D N ^ 
(thus N < N). Then from © 

C k - n D k = for any n < N 

k>N;k^0 

or, setting fc = TV — 1 + j, 

Y C j+n D N - 1+j = for n > (32) 

j>i,jyi-jv 

It is here that we use the genericity condition on C. In fact we will show that 
([52|) implies D = if condition (g) is satisfied. To see this define D € h(N) as 
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Dj = D N _ 1+j if 3 > l,j ^ 1 - N and, if 1 - N > 1, Di_;v = 0. Then by © 
£) is orthogonal in Z 2 (N) to all T n C, n > 0. By the genericity condition (g) then 
< D, e\ >= Dm = 0, which is a contradiction. Thus D = 0. 

Since J' is analytic in sq for small enough sq, and compact by the same 
simple arguments as in Lemma ^ it follows that (I— J')~~ x exists and is analytic 
in so at sq = 0. 

(ii) This part is an immediate calculation. 

(hi) Note first that f £ h{^ \ {0}), because 



< 



l - VTT^ 



2s 




E 



Ck—n 



Sq + kiU 



1/2 



<iwie- 



•so 



kuj\ 2 



< CO 



Also, / being an integral with respect to a discrete measure of an l\ analytic 
function, depends analytically on small sq. 

The rest of the proof of (hi) closely follows that of part (i), using the following 
result. 

Lemma 11 For s — we have Cjqjdj ^ 0. 

Proof. Assume the contrary was true. At sq = 0, with D° = D n \ So= o and 
q° = 3nL=o, relation @, using @, gives 

D° 



9 " fe#0 



(33) 



Multiplying with D° and summing over we would get 

^(-1 + VI+^)K| 2 = - J2 C k - n D° k Dj-^C-nDj (34) 

and since we assumed ^ n C n D^ = then, as in the proof of Lemma [lO] (i), it 
follows that D° = for all n < N = -W 1 . This gives, using (||), that 



2 Cfe-n£>fe + 2C_„ = (35) 

fc>Af o; A:^0 

Denote by D 1 £ I2 the sequence D\ = if k ^ and Dp = 2. As in the 
proof of Lemma [Ic] (i), using the genericity condition (g), we get D 1 — 0, an 
obvious contradiction. 

This concludes the proof of Proposition [| for generic 77 the solution y of ( |l7|) 
has, for lu^ 1 ^ N, analytic components y n when p = 0. 

Square root singularities. We now study the behavior at the square root 
singularities of the coefficients of the equation of y. 

Let fco be the unique integer such that for some s r £ [0, lo) we have 1 + ,s r + 
fcow = (then s r is a branch point in the coefficient q). 

The following Proposition describes the analytic structure of y(p) near the 
imaginary axis. 
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Proposition 12 We have the decomposition y n = u n + {\Jsq — s r )v n where u n 
and v n are analytic in sq in a complex neighborhood of the segment [0, u>). 



Proof. The substitution y n = u n + (-\/io — s r )v n , and 

Uk = qku k ; Vfc = q k v k (k ^ k ) and U ko = — — ; V ko 



SO + k LL> s + koto 

leads to the following system of equations for U n and V n 



q n l U n = riV — - " y^Ck-nUk - C ko - n (so ~ s r )V ko (n 7^ fc ) 

3n l V n = — } j C k ~ n V k - Cfc _„(s - s r )V ko - C k „- n U ko (n ^ k Q ) 



k 



Ck- 



k u 



(s + k u)U ko = i V] 

^ s + few 

A: 

(s + fcow)T4 =~J2 C k-kaVk (36) 



We now let Q k(> = s + k uj and, for n ^ k n , Q n = q~ x = -1 + vT+so + kuj. 
We use again the Fredholm alternative and, as in the previous proofs, we need 
only to show the absence of a solution of the homogeneous equation at sq — s r . 
We thus multiply the homogeneous equations associated to (|36|) in the following 
manner: the equation for Uj by Uj and the equation for Vj by Vj, then sum over 
all j. As in the previous proofs, from the reality of the r.h.s. and then from the 
genericity condition (g) U = 0. Then, similarly, V = 0. The rest is immediate. 



5.2. The resonant case: lo^ 1 = M G N. In this case when so = there are poles 
in the coefficients of ( p3| ) when n + j = and branch points when n + j = —M. 
The proof is a combination of the two regularization techniques used in the 
previous case. 

Proposition 13 We can set y(so) = A(sq) + B(so)y/s~o with A and B analytic 
in a complex neighborhood of the segment [0,o;). 



Proof. Special care is only needed near so = 0. The system (E6|)-(28) now reads 



d n — —C- n (l + VI + So) — ^ C k - n q k d k — C-M-n —d-M 

— so — 1 

k${0,-M} u 



v n — fn — ^ Ck- n qk v k ~ C-M-n ~ V -M (37) 

so — 1 

ki{a-M} ' 
We take d n = ot n + ^/soPn and v n — j n + ^fs~o8 n . The system becomes 
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a n = -C- n (l + VT+ S ) - Ck-nQkCtk- C-M-n ~ (OL-M + S f3-M ) 

So — 1 

H{0-M} 

Pn = — ^ Ck-nQkPk — C-M-n ~{ a -AI + P-Al) (38) 



s - 1 



7n = /n — ^ Ck-n<lklk - C-M-n ~ (7- Af + So<5-j\f) 

. . . _ So - 1 



so 

fc£{0,-M} 



0~n = - \^ Ck-n<lA - C-M-n ~{5-M + 7-Af) (39) 

— So — 1 

fc£{0,-A/} U 



The system 138) is of the form 



<A s{so) M , m 



where a,f3,F\,F% are in ^. We prove that the homogeneous equation has no 
nontrivial solutions: 

Lemma 14 (I - 5(0)) (j^j = implies r^j = 0. 

Proof. Let Q„ = q n , A n = q n a n , B n = q n f3 n for n ^ 0,—M and Q-m = -1, 
= — a_M and £?_m = —0-M- The system (g8j) becomes 

Qn — — ^ ' Cfc— n Af~ 

k^Q 

Qn Y Bn = — 2^ Ck-nBk — C-M-nA-M (40) 
k^O 

As in the proofs in Case I, multiplying the first equation by A ni summing over 
n we first get from the reality of the r.h.s. that A n = for n < —M and then by 
the condition (g) we get that A = 0. The conclusion B = now follows in the 
same way. 

End of proof of Proposition [7^. The operator 5 is compact on l 2 © h and 5 and 
(Pi, F%) are analytic in a complex neighborhood of 0. We saw in Lemma [l4|that 
the kernel of I — 5(0) is trivial and by the analytic Fredholm alternative it follows 
that (I — 5(0)) _1 exists and is analytic in a small neighborhood of so = 0. Hence 
(a, 0) are analytic. Similarly, j,S are analytic in the same region. 



5.3. Proof of Theorem |J. Combining the above results we have the following 
conclusion: 

Proposition 15 If condition (g) is fulfilled, then y(jp) is analytic in a neighbor- 
hood of iR \ {is r + iwZ}. For any j € Z, in a neighborhood of p — is r + ijoj 
(s r € M.) y has the form y{p) = Aj(p) + Bj(p)y/—ip — s r — ijuj where Aj and 
Bj are analytic. In particular, y is Lipschitz continuous of exponent 1/2 in the 
closed right half plane. Thus Y(t) = 0(t~ 3 / 2 ) for large t. 
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Proof. All but the last claim has already been shown. The last statement is a 
standard Tauberian theorem (note that C^ 1 is the Fourier transform along the 
imaginary line). 

Proposition 16 We have 9{t) — > as t — > oo. 
Proof. We can write @ (with I(t) = 1) as 



Y = r](9 + M*Y) (41) 

It is easy to check, in view of the fact that M and Y are 0(t -3 / 2 ), that M*Y —* 0. 
Furthermore 1 + 2i J Q Y(s)ds is convergent as t — > oo. Thus 8(t) — > const as 
t — > oo. Since now the l.h.s. of (Q) converges to zero and 77 does not, the equality 
( [HI) is only consistent if — > 0. 

This completes the proof of Theorem [| for the case tpo = Ub = e - '^'. 

The general case follows by noting that the inhomogeneous term does not 
affect the main argument, using the Frcdholm alternative. Hence we will still 
have \0(t)\ — > but the rate of decay may be different. 



6. A nongeneric example 

Let 77 be given by (pi]), for which 

C n = -r\ n forn>l , C n = C_„ (42) 
As in § 5 set —ipo = so and let q n be given by (Pj). Denote 

a n = a n («o) = -— = - (Vl + s +nuj - l) (43) 

For re (0, 1), ui > 1, ^ N such that (1 — r) 2 < a; — 1, let s r and s p be the 
unique numbers in (0, ui) so that 1 + s r G luZ and 1 + a-i(s p ) = 0. We choose 
r, w such that s r ^ s p . 



6.1. The homogeneous equation. 

Lemma 17 Lei so,o be a point in (0, s r )U (s r , w). Consider so m a small enough 
neighborhood o/so,o- 27ie linear operator J — J {so) of (Ujj depends analytically 
on so, and is compact on ?2- For so 7^ s p , (I — J (so))~^exists and is analytic. 

Lemma 18 Denote for short Jo = J(s r ). 

There exists a value A = A s € (0, 1) such that 

dim Ker (I - Jo) = 1 (44) 

Denote by A the diagonal (unbounded) operator (Az) n — a n z n in I2', A^ 1 is 
bounded. 

Lemma 19 For A = X s as in Lemma we have 

Ker (I - Jo) = A[Ker(I - J*)} (45) 
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6.2. Proof of Lemma |7J. The operator J is compact by Lemma |^. To show 
that / — J is invertible we prove this for any points sq G (0, w), sq ^ s p , s r ; by 
the analytic Fredholm theorem it will follow that / — J is invertible in a small 
enough neighborhood of any such point, thus proving the Lemma. 

Let so € (0,o->), so ^ s p ,s r . As in Remark (?] in sectio n [s| , the substitution 
y n = a n z n (n £ Z) transforms the homogeneous equation (|19|) to 

oo 

a n z n = ^2 A J (z n+ j + z n -j) ■, n eZ (46) 
i=i 

Note that ^a n < for n < — 1 for so £ [w — 1,oj) and 3a„ < for n < for 
so G (0,w — 1). We will discuss only the first case, sq > u> — 1, since the second 
one is completely analogous. 

As in the proof of Lemma ||, it follows that 

z n = for n < -1 (47) 

Then equations ( ^6|) for n < — 1 become 

oo 

£ A fe z fc _2 = (48) 



k=i 



For n. = — 1 (Eq) gives 



(o_i + l)z-i = (49) 
and for n > 0, using (f48|), we get 

n+l 

(1 + a„)z„ = ^(A J - *~ j >n-i > n > -1 (50) 

3=1 

Since so ^ s p , ( (49|) gives Z-i — 0, and it follows by induction, from (fiO|), 
that z n — for all n. By the Fredholm alternative theorem then / — J'(so) is 
invertible. 



6.3. Proof of Lemma^T^. In what follows so = s r . 

6.3.1. An auxiliary lemma. We show that if z 6 I2 then equation (|4§| ) is redun- 
dant. 



Lemma 20 //z is an ^2 solution of ( ]5(\ ) with z n — /or n < — 1 then z satisfies 
&■ 

Proof. Let z £ I2 be a solution of (|50|). Then 

n 

z [«+i] = J- \ k z k _ 2 (51) 
fe=i 



18 



O. Costin, R. D. Costin, J. L. Lebowitz, A. Rokhlenko 



is the truncation of a convergent series, since there is a constant M with \z n \ < M 
for all n. Note that 

n+l 

(1 + a n )z n = A^n-i - A-"- 2 Z[" +1 l 

3=1 

hence 

n+l 

Z [n+1\ = X n+2 \j Zn -i - A"+ 2 (l + a n )z n 



so that 

Z [n+1] 



< A n+2_^_ + x n+2 M(l+ const Vn) — >0 asmoo (52) 
1 — A 



Since ( (Slf ) are truncations of the series in the LHS of (p8[), then (p2| ) implies 



6.3.2. Behavior of the general solution of §5(\). A direct calculation shows that 
the sequence z n satisfying the infinite order recurrence ( |50| ) and the initial con- 
dition z_i = 1 satisfies, in fact, the three step recurrence 



(1 + a n+ i)z n+1 + (1 + a„_i)z„_i 
with the initial condition 



[A(l +a n ) + A + a„A" 1 ]z„ (n > 0) (53) 



z-i 



1. 



^0 



A -A' 1 
1 + a 



Denote 



then (p3t) becomes 



A-A-\ , 

Zn = -— V n -1 



V n + V n - 2 = 



A + 



A 2 + a n 



A(l 



n > 1 



(54) 



(55) 



(56) 



We are looking for l 2 solutions. Standard WKB (see e.g. Q) would imply 

there are solutions of the discrete equation ( plf ) behaving like A _n e _ v / ™/" and 

like A™eV / "/". We will prove this in our context and find special values of A for 
which the solution decaying for large n satisfies the initial condition. We will 
show that this solution is obtained by taking 



Vn-2 — 9n-lVn-l 



(57) 



in (56) and iterating: 
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i.e., go is given by the continued fraction: 

1 



.go = G\ - 



which needs to match the initial condition (see ( p4j 

1 

A + A- 1 + (l + a )" 1 (A - A- 1 ) 



5o = 9o (A) = ~ . , t~ ~~ -— r- (59) 



Lemma 21 (%) .Lei A G (0,1). The recurrence ( |5q j has a solution such that 
g n — > A -1 as n —>■ oo. 

(ii) go is meromorphic in A on [0, 1) and has poles. 

(Hi) There exists X s € (0, 1) swc/i i/iai <7o(A s ) satisfies $5(\). 

('ivj Lei A = A s . To i/ie solution of (i) there corresponds a solution of the 
recurrence $5(\ ) such that Vn^ ~ As^ '™' as n — > oo. TTie corresponding solution 
z^l o/ ( pZ^ satisfies z n — > as n — > oo. 

(V) Lei A = A s . There exists a solution of toztft. i/ie asymptotic behavior 

y[l] ^ ^-n+o(n) 

Thus, for A = A s , i/iere exists a unique (up to a multiplicative constant) 
"small" solution of $h\), with the behavior Vn^ ~ A™ + °^ /or /arge n, w/ii/e 
£/ie general solution behaves like V n ~ A s n+ . vis a consequence, a similar 



statement holds for the recurrence 

Remark. The proof of (hi) can be refined to show that, in fact, there is a 
countable set of points A s for which go satisfies the initial condition, and that 
these values accumulate to 1. 

Proof, (i) With the substitution 

9n = G n+ i - A + 5 n (60) 



the recurrence (pi) becomes 



G n +2 — A + 6. 



71+1 



(SS) n , n>0 (61) 



For no > and e small, positive, define A„ = a„ 0+ 2 (2 + a„ +2) _1 — £• Let 
Af no be a small neighborhood of the interval I no = [0, A„ ]. Consider the Ba- 
nach space B no of sequences {S n } n>n with S n — S n (X) analytic on Af no and 
continuous up to the boundary, with the norm = sup Jl>Jlo sup^g^ |5„(A)|. 
Direct estimates show that the operator S defined by (^lj) takes the ball of radius 
p no — 2/(2 + a„ 0+ 2) + e' in B no into itself (if e,e' and M no are small enough), 
and is a contraction in this ball. Therefore the equation 5 — S(S) has a unique 
solution in B no , of norm less than p„ . Then |o~ ra (A)| < const(n + 2) -1 / 2 for all 
A € 7 n and all n > 0. Since the sequence A„ increases to 1, (i) follows, 
(ii) Step I: all g n are meromorphic on [0, 1). 

Since S n is analytic on /„, then from (pw, g n is analytic on /„ \ {0}, having a 
pole at A = 0: g„ ~ A _1 a n+ i(l + a n+ i)~^(A — > 0). Iterating (Eg) it follows that 
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<7n_i, g n ~2, ■■■,go are meromorphic on I n . Since the intervals /„ increase toward 
[0, 1) it follows that go,9i, ■■■9n--- are meromorphic on [0, 1). 

Step II: there exists n\ and Ao € (0, 1) such that 3m (Ao) < 0. 

Define e„ = (1 + a n )~ ; we have (see J43|)) 

e« ~ r(n uj)~~ 1/2 , n —> oo (62) 

Let no be large and denote Ao = 1 — e rao . Let No be large enough so that Ao 
is in the domain of analyticity of gN - Iterating fl58| ) starting from Nq (and 
decreasing indices) we get the value <M (Ao). If for some n £ {rig, uq + 1, Nq} 
we get g n (Xo) < 0, Step II is proved. Then assume that g„ (Ao) > 0. 
Consider the recurrence 

f n _i = G„ (A ) - — for n < n , r n = 5„ (A ) (63) 



2 

no ' 



where, in fact, G no (Ao) — 2 — e 

The recurrence ( p3| ) can be solved explicitly (it is a discrete Riccati equation 
and a substitution f n — l = x n -i/x n transforms it into a linear recurrence with 
constant coefficients). It has the solution 

cos((n- n )(p + 6) 
■ - (64) 



cos (in + 1 — n )</> + 0) 

where cos0 = 1 — e„ /2, sin^i > 0, and tan# = (cos0 — X)/sm(j) so that 9 ~ 

5 ~ i e «o ( e "o — * o)- 

We assume, to get a contradiction, that ^n(Ao) > for all n = 0, 1, ...,m. 
Then 

g n {Xo) < f n for n <n (65) 

which follows immediately by induction using (^8|), (|63|), noting that G„ is in- 
creasing in n. 

Note that there is an n\ £ {1, 2, no — 1} so that 

f n > for n G {n x + 1, ...,n } and f ni < (66) 

Indeed (from (|62])) when n decreases from no the numerator and denominator 
in ( |64| ) increase up to 1, then decrease, until the numerator becomes negative, 
when n equals n\ — no — k\ where k\ is the integer with k\ — 1 < (w/2 + 9)/<f> < 
ki. Since ~ e„ n (e„ — > 0) then k% ~ (37r)/(4e„ ), and, using (|62|), clearly 
fci £ {1, rtp — 1} (if no is sufficiently large). 

Then ( p5| ) and ( |6q ) contradict the assumption that <7 ni (Ao) > 0, and Step II 
is proved. 

Step III. 

The function g ni is meromorphic on [0, 1), with g ni (0+) = +oo. There is a 
smallest value of A in (0, Ao) where g ni changes sign: this is either a zero, or a 
pole. 

Assume it was a pole. Let p € (0, Ao) be the first pole of g ni . Then g ni 
is positive and analytic on (0,p), and g ni (p~ ) = +oo, g ni (p+) = — oo. Since 
g n+ i = l/(G n+ i - g n ) (see @) then g ni +i(p-) = 0-, hence g„ 1+ i changes 
sign in (0,p). But g ni +i has no zero in (0,p) (otherwise at that zero g ni would 
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have had a pole, from (|5q)). Then g ni +i has a pole, with a change of sign, from 
+ to — , in (0,p). Now the argument can be repeated. It follows that for any 
k > 0, g ni +k has a pole in (0,p), which contradicts the fact that the domain of 
analyticity of g n increases to (0,1) as n — > oo. 

Therefore, the first change of sign of g ni is at a zero. Let (i be the smallest 
value in (0, 1) such that g ni (Ci~ ) = 0+, g ni (&+) = 0—. Then from (58) we have 
<7m-i(Ci — ) = ~°° an( i 9W-1 changes sign in (0, £i)' Now the argument can be 
repeated. It follows that go has a pole at a point £ ni with <7o(Cni — ) — — oo. 

(hi) Since flp(A) takes all the values when A S (0, £ ni ) there exists A = X s e 
(0, 1) such that (||) holds. 

(iv) For A = A s , since the solution of (i) satisfies g n (X) = A -1 + 0(n -1 / 2 ) we 
have from @,with the notation = V(X a ), that V,[ s] = flLo fffe(A a )~V [sl = 
0(Ar o(n) ) and thus V^-V^ = 0(A™ +o(n) ); then from © zjf 1 = 0(A s l+o(n) ). 

(v) The substitution (variation of constants) V n = V-h v n brings the recur- 
rence (|56l ) to a first order one: with the notation A n = v n — v n -i we have 

A n = V^[_ 2 /^n An— i and the rest of the argument consists of straightforward 
estimates. 



6.3.3. Proof of Lemma |Z§ 



Proof. Lemma |2l|(v) shows that equation (53) has a unique (up to a multiplica- 



tive constant) small solution, ~ A" + ° < "' ) (n — ► oo), while the general solution 
behaves like z n ~ \ s n+ °( n \ Since y„ ~ \fnz n the uniqueness of the I2 solution 
is proven. 



6.3.4- Examples of solutions. We will show next how concrete values A s satisfy- 
ing Lemma |2l| (hi) are relatively straightforwardly, and rigorously, found. One 
method is as follows. Note that the minimum/maximum of the function a — b/x 
where x varies in an interval not containing zero is achieved at the endpoints. 
We thus take the recurrence (|5q) with initial conditions g no = A -1 ± ^= and 
compute go from these. The actual graph will be between these two, unless the 
condition mentioned is violated in between no and 0. This graph is to be inter- 
sected with the graph of the initial condition ( p9| ) . 

We take for instance w = 1.1, r = 0.45, s p = 0.11, n = 10, for which 
the rigorous control is not too involved. The two graphs are very close to each 
other (within about 3.10 -6 for A S (0.3,0.4)) and cannot be distinguished from 
each-other in Figure 1. A first intersection is seen at A ~ 0.327; see Figure 2. 



6.4- Proof of Lemma [/I]. Denote B = (I — J^o)A; we have B = A — S. Hence 
B* = ~A- S. Then Ker(B) = Ker(S*) (since Az = Sz implies @, so Az = Az, 
and similarly, Az = z implies Az = Az). So Ker[(l - Jo) A] = KerL4(l - Jo*)] 
so that (since A is one-to-one) A _1 Ker(l — Jo) = Ker(l — Jo*), which proves 
the Lemma. 
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Fig. 1. Graph of go given by ( pq ) (discontinuous graph) and by ( |59| ) in a region near A = 1, 
as functions of A. 
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Fig. 2. Graphs of go (steeper graph) and of the initial condition for go ( 
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6.5. Discussion of the singularities of solutions of Let A = A s . We have 

that I — J is invertible for 5Rpo > 0, and is not invertible at po = is p (Lemma 
|l8| ). By the analytic Fredholm theorem (see e.g. |p0| ) (I — J)^ 1 is meromorphic 
on a small neighborhood of is p , therefore there exist m > 1 and operators 
S m , Si,R(p ) so that: 

(I-JV 1 t 1 —^S m + ... + — Si + R( Pa ) (67) 

(pa-is p ) Po-isp 

where R(po) is analytic at is p , and S m ^ (since I — 3§ is not invertible). 
Multiplying (|67|) by / — J to the left, respectively to the right, and writing 
J = J® + (po — «s p )j7i(po) (where Ji(po) is analytic at is p ) we get that 

Rl(Po) = 7 ym {I ~ Jo) S m + O ({po ~ is p y m+1 ) 

(Po-is p ) V / 

R2(Po) = 7 ~. ym Sm i 1 ~ Jo) + ({p - is p y m+1 ) 

(PO -ISp) V / 

where R± t 2 are analytic at po = is p . By the uniqueness of the series of the 
analytic functions (Banach space valued) Ri,2 we must then have 

(I-Jo)S m = = S m (I-Jo) (68) 

The first equality in ( |68| ) implies Ran(S' m ) C Ker (/ — J ) = \/{yKer} and 
since Sm 

^ then Ran(S' m ) = \J{yKer} therefore S m y — (y, u)yKer for some u £ 
l 2 \ {0}. The second equality in ( p8| ) means u £ Ran (/ — J ) = Kcr (/ — J~ *). 

By Lemma [l^ then (up to a multiplicative constant) u = A~ l yxer = ^Ker 
where z^ er satisfies (p6|), hence (|53|) , (|54|) . The solution y = (I — J) 1 f o f (E~ 



is then singular at po = if c =< f,ZKer 0- For the example of §6.3.4 
this latter condition can be checked by explicit calculation of the truncations to 
10 terms and estimation of the remainder based on the contractivity bounds in 
the previous section. The result is c = — 1.953 ± 0.001. Thus the inhomogeneous 
equation has poles. 

Lemma 22 Let Y(t) be analytic on [0, oo), with limt—nx, Y(t) = 0. 
Let s6l Then 

lima/ e- {a+is)t Y{t)dt = (69) 
<40 J 



Corollary 23 Let Y(t) be as in Lemma^Q. Let y{p) = / °° e~ pt Y(t) dt. Assume 
that y(p) is analytic on £R+, except for a set of isolated points. Then y(p) does 
not have poles on 

Proof. L. We first show (|9|) for s = 0. 

Separating the positive and negative parts of $tY(t), ^sY(t) write Y(t) = 
- Y^(t) + iY^(t) - iyW(i) with FW(t) nonnegative, continuous, non- 
analytic only on a discrete set, where the left and right derivatives exist, with 
limt^oo Y^(t) = 0. It is enough to show ([39]) for each Y^. Let then Y be one 
of the yH's. Denote H(t) = sup T > t F(r). The function H on [0, oo) has the 
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same properties as Y and in addition is decreasing. Then H' exists a.e. and 
H 1 e L x [0, oo) since /°° \H'(t)\ dr = - lim*-oo /„' H'{t) dr = lim*-oo -H(t) + 
H(0) = H(0) 
Then 

poo pOO pOO 7 

a / e- at F(t) dt < a / e- at H{t) dt = - — {er at ) H{t) dt 
Jo Jo Jo 

poo 

= H{0)+ / e~ at H'(t)dt 
Jo 

therefore 

/>oo />oo 

lima/ e- at Y(t) dt < H(0) + lim / e~ at H'(t) dt = 

"10 Jo aj.0 ,/ 

which proves the Lemma in this case. 

77. Let now s £ R arbitrary. Then (^) follows from the result for s = 
applied to the function Y(t) = er lst Y(t). 

Proof of Theorem |2] 

In conclusion Y(t) cannot tend to zero as t — > oo and complete ionization 
fails. 
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Appendix A. Proof of Lemma g 

(i) Consider Lj oc [0, A] endowed with the norm ||F||„ := J \F(s)\e- vs ds, where 
v > 0. If / is continuous and F,G £ Lj oc [0,A], a straightforward calculation 
shows that 

||/ J F|U<||F|Usup|/| (Al) 

[0,A] 

\\F*G\\ V < \\F\\ U \\G\\ V (A2) 
\\F\\ V -^0 as v -> oo (A3) 

where the last relation follows from the Ricmann-Lcbcsguc lemma. 
The integral equation (|^) can be written as 



Y = r] + JY where JF := ri{2i + M) * F (A4) 

Since M is locally in L 1 and bounded for large x it is clear that for large enough 
B2, (§) is contractive if v > B2, for any A. 

(ii) This is an immediate consequence of Lemma [| and of elementary prop- 
erties of the Laplace transform. 
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(iii) We have in 



2* [ , -r,r f u 2 e - l( x- 4 a)(l+„-) 

/ Axe px / - 

alO 7T J J 1 + li 2 



CM = lim — / dxe~ px / rT - 5 du (A5) 

du (A6) 



W-co (1 + u 2 )(p + i(l + u 2 )) 



For 3?(p) > we push the integration contour through the upper half plane. 
At the two poles in the upper half plane u 2 + 1 equals and ip respectively, so 
that 



i 



oo „2 

du 



W-oo (l + u 2 )(p + i(l + u 2 )) 



(-1) /ds u 2 , i^M = _l + I^ (A7) 



7r \(2i)(p) J s (ip)(2iuo) J s J p p 
where uq is the root of p + i(l + u 2 ) = in the upper half plane. Thus 

CM = - i - + i ^^- (A8) 
P P 

with the branch satisfying \fl — ip — > 1 as p — ► in BL As p varies in H, 
1 — ip belongs to the lower half plane —ill and then ^/l — ip varies in the fourth 
quadrant. 

For Sft(p) > 0, uj > we have 



C(e^M) = *- + ^- ip ™ (A9) 

/ p =F icj p^f itu 



(with yl — ip — a; = —iy uj — 1 + ip if a; > 1) 
and relation (Il5|) follows. 



Appendix B. Discussion of the genericity condition (g) 

A thorough analysis of the properties of the shift operator is provided by the 
treatise pil. We provide here an independent discussion, meant to give an insight 
on the interesting analytic properties involved in this condition. 

Let C = (Co, Ci, C n , ...) € ^(N) and the operator T defined as before 
by TC = (Ci, C2, ■■■)• We want to see for which such vectors, the system of 
equations 



(z,T'C) = 0, J = 0,1,... (Bl) 

has nontrivial solutions z in /2- We can associate to z and C analytic functions 
in the unit disk, Z(x) and C(x) by 
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C(x) = J2C k x k Z(x) = J2 



z k x 



(B2) 



k=0 



k=0 



These fun ctions, extend to L 2 functions on the unit circle. The system of equa- 
tions (Bl) can be written as 



z C(x) + Zxx- 1 {C{x)C{!S}) + ... 



71—1 h 

X 



- n C(x)-x- n J2lJ c{k) ( ) 



k=0 



= (B3) 



Using Cauchy's formula, we can the difference in square brackets in(B3) as 



C(s) 



2tti 7j s]=1 s n (s - x) 



-ds 



(B4) 



and thus (|Bl|) becomes 



C{s)Z{l/s) 



ds = 



(B5) 



| s |=i 



The functions C for which this equation has nontrivial solutions Z relate to the 
Beurling inner functions (29) and are very "rare" . 

Examples, (i) Let |A|< 1 and C n = A™, i.e. C(x) = (1 - Xx)' 1 . This is 
related to the function (jO). Taking advantage of the analyticity of Z outside 
the unit circle, we can push the contour of inte grat ion towards s — oo, collecting 
the residue at x = A -1 ; we see that equation (B5) holds iff Z(X) = 0, i.e., for a 
space of analytic functions of codimcnsion one. 

(ii) Let X n = 1/n. Then C(x) = ln(l — x), and by taking s — l/t in ( |B5| ) we 

get 



1 / Z(t)]n(t-1) 1 / ln(t)Z(t) , , 

- & )' y ,' dt - - i , y ' \\ dt = B6 
xj\t\=i (t-x-i)t xj w=l t{t-x-^) v ' 

By making a cut on [1, oo) for the log we see that the integrand in the first integral 
is analytic in the unit circle and thus the integral vanishes. We decompose the 
second integral by partial fractions and we get 

r mm dt = (B7) 

T\t\=i t J\ t \=i (*-y) 

where y = x^ 1 . The first integral is a constant, C. By pushing the contour 
of integration inwards, we see that the second integral extends analytically for 
small y ^ 0. For such y we thus have 
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Hmz(t)-z(y)) dt + z{y) r Mt) dt = _ c m 



1*1=1 (t-y) J\t\=i(t-y) 

Now the contour of integration can be pushed to the sides of the interval [0, 1] 
collecting the difference between the branches of the log. We get 

^-^WW 1 — <** = (B9) 



t-y Jo t-y 



Thus <fi(y) + Z(y) ln(— y) — C with <fi and Z analytic in the unit circle thus 
ln(— y) is analytic unless Z = 0. This shows C n = 1/n is generic. 
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